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Building on previous work on the critical behavior in gravitational collapse of the self-gravitating 
SU( 2) (T-field and using high precision numerical methods we uncover a fine structure hidden in a 
narrow window of parameter space. We argue that this numerical finding has a natural explanation 
within a dynamical system framework of critical collapse. 


Introduction 

Over the past few years the Einstein- SU{2) — a model has attracted a great deal of attention 111111110 . This 
model is interesting because its rich phenomenology is sensitive to the value of a dimensionless parameter rj which leads to 
various bifurcation phenomena. The most interesting bifurcation was found by Lechner et al. 0] who showed that the critical 
behavior in gravitational collapse changes character from continuous to discrete self-similarity when the coupling constant 77 
increases above a critical value r] c . This phenomenon was interpreted in terms of dynamical systems theory as the homoclinic 
ioop bifurcation where the two critical solutions, continously self-similar (CSS) and discretely self-similar (DSS), merge in 
phase space. Since the echoing period A of the DSS solution diverges as 77 —> r] c , the numerical analysis of this bifurcation 
is extremely difficult and for this reason some of the aspects of critical behavior near the bifurcation point, in particular the 
black hole mass scaling law, were left open in Q. 

Below, using high precision numerical methods, we confirm the main findings of 0. In addition, we find that just above the 
bifurcation point the marginally supercritical side of the transition between dispersion and black holes exhibits a fine structure 
which is due to the competition between two coexisting critical solutions, the DSS one and the CSS one. The description of 
rhis phenomenon and its interpretation is the main purpose of this paper. The rest of the paper is organized as follows. For 
readers’ convenience, in section 2 we first briefly repeat the basic setting of the model and then we summarize what is known 
about it. In section 3 we present numerical results and finally, in section 4, we interpret them. 


The model 


The spherically symmetric Einstein-5'/7(2) — a system is parametrized by three functions: the metric coefficients A(t,r), 
6(t,r) and the cr-field F(t,r ), which satisfy the following system of equations 




U g = —e A dt(e A A~ 1 dt) + -=■ d r (r 2 e~ d A d r ) 


(1) 


d t A = -277 rA(d t F)(d r F), (2) 

d r 5 = -rjr({d r F) 2 + A- 2 e 2 \d t F) 2 ), (3) 

d r A = l^-A- r] r(A(d r F) 2 + A- 1 e 25 (d t F) 2 + 2 V (4) 

r \ r z J 

where 77 is a dimensionless coupling constant. For 77 = 0 this system reduces to the a model in Minkowski spacetime. The 
initial value problem for this system was studied by Bizon et al. '9] for 77 = 0 and by Husa et al. jj for 77 > 0. In these studies 
an important role is played by self-similar solutions. A countable family of continuously self-similar (CSS) solutions, herefater 
denoted by CSS n (n = 0,1,...), was shown to exist for 0 < 77 < 0.5 in illS These solutions are regular within the 
past light cone of the singularity, however they have a spacelike hypersurface of marginally trapped surfaces, i.e. an apparent 
horizon outside the past light cone if 77 > r ] n , where ij n is an increasing sequence (770 = 0.08, 771 = 0.152, etc.). Linear stability 
analysis shows that the ’’ground state” CSSq is stable while the excitations CSS n have exactly n unstable modes. 

Besides the CSS solutions, the system (1-4) has also a discretely self-similar (DSS) solution for 77 > 77 c ~ 0.17. This solution 
was constructed by Lechner ( 3 ] via a pseudospectral method following the lines of Gundlach |j. 

Next, we summarize what is known about the critical behavior in gravitational collapse in this model. The first numerical 
studies of this problem, reported in Q, focused on relatively large coupling constants 77 > 0.2. In this range a ’’clean” type 
II critical DSS behavior was observed, however the attempts to resolve critical evolutions for lower values of 77 encountered 
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numerical difficulties and for 0.18 < 77 < 0.2 only an approximate DSS behavior was observed. Furthermore the echoing 
period A was found to increase sharply as the coupling constant decreases from 0.5 to 0.18. The critical behavior for smaller 
couplings 0.1 < 77 < 0.2 was studied in Q (still smaller couplings are less interesting because then the model admits naked 
singularities). In the range 0.1 < 77 < 0.14 a ’’clean” CSS critical behavior was observed, thus it became clear that somewhere 
in the interval 0.14 < 77 < 0.2 there must be a transition between CSS and DSS critical solutions. The detailed studies of 
this transition 3 led to a conjecture that there exist a critical value of the coupling constant r) c « 0.17 for which the system 
exhibits the homoclinic loop bifurcation, i.e. the CSS saddle merges with the DSS limit cycle in the phase space. These results 
left open the question which of the two solutions in the transition region acts as the critical solution at the threshold of black 
hole formation. In particular, near the bifurcation point the black-hole mass scaling could not be properly resolved. 


Numerical results 

We have solved equations (l)-(4) for marginally critical initial data fine tuned to the DSS solution for coupling constants 
close to the critical value rj c = 0.17. Since the echoing period A increases sharply as the coupling constant tends to its critical 
value from above, it becomes more and more difficult to follow the evolution over a large number of DSS cycles 1 . We used the 
fully constrained implicit evolution scheme based on the Newton-Raphson iteration. In order to resolve the singular behavior 
near the origin we used the grid which is uniform in ln(r). To get several cycles of the DSS attractor near the bifurcation 
point we had to fine tune parameters of initial data with precision of 70 digits - this was achieved with the help of the arprec 
library pj| 

Actually, it was not our aim to determine A with high precision, but rather to show that there exists a A to which the 
evolution converges. To this end, we determine A as a function of time (cycles) as the evolution approaches the limit cycle i.e. 
the DSS solution. For a marginally critical solution we plot the function F versus ln(r) for some late time t.\ and superimpose 
the profile of the first echo at time £2 shifted by ln(?-) —> ln(r) + A r . The time t% and the radial echoing period A r are 
chosen to minimize the discrepancy between two profiles. We also define the temporal echoing period A f by the formula 
t ‘2 = t*( 1 — e -At ) + e -At fi, where t* is the accumulation time. Repeating this calculation for a sequence of pairs (t n , t n+ 1 ), 
we get a sequence of values A r and A t . Of course, if the evolution converges to the DSS solution, both A r and A t should 
converge to the same constant. In Fig. 1 we show the convergence of A r and A t during a critical evolution for the coupling 
constant 77 = 0.1725. Note that the curve levels off, thus signaling the closeness of the evolution to the limit cycle. As the 
coupling constant decreases, A grows and the approach to the DSS solution becomes slower. 



0 2 4 6 8 10 


N 


FIG. 1: For rj = 0.1725 we show A r and At as the functions of the cycle number N. Fitting the curve A + ce to the data we obtain 
A ss 1.803. 


1 By an elementary dimensional analysis the number of cycles scales as n ~ — (l/AA)ln|p — p*|, where A is the echoing period and A is the 
eigenvalue of the growing mode of the DSS solution. 
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FIG. 2: Fitting the 
and the slope a = — 
H] Acss(??c) ~ 5.14. 


echoing period A (determined as in Fig. 1) to the analytic prediction A 
0.389 which is in very good agreement with the analytic prediction a = — 


= a In \q — rjc \ + const 
2 /A css and the linear 


we get rjc = 0.1701 
perturbation result 


Let p* be a critical parameter value which separates dispersion from black holes (this value can be found by standard 
bisection). In agreement with [I] we find that for ?y > 0.17 the solution corresponding to p* is DSS, in particular for p slightly 
below p* we observe DSS subcriticality, i.e. the solution approaches the DSS solution and then disperses. Looking at the 
maximum value of the spatial derivative of the scalar field at the origin as the function of p, we find a typical subcritical 
scaling law (see Fig. 3) 

max \d r F(t, 0)| ~ (p* —p) _7DSS . (5) 



FIG. 3: The subcritical scaling (5) for the coupling constant 77 = 0.19. The slope of the linear fit is approximately equal to —0.109. The 
wiggles, which are imprints of discrete self-similarity, have the period « 4.8. 

For p > p* black holes are formed, however this happens in a rather unusual manner. This is shown in Fig. 4 where the 
metric function A is seen to develop two minima very close to zero which signals an almost simultaneous formation of a small 
and a large black hole. 
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FIG. 4: The series of snapshots of the metric function A(t, r) from the evolution of marginally supercritical initial data for the coupling 
constant rj = 0.19. At the last frame one can see the formation of two apparent horizons. 


Let us denote their apparent horizon radii by rj n and r out , respectively. We find that the outer radius exhibits the standard 
DSS supercritical scaling (see Fig. 5a) 


router ~ (P ~ P*V DSS , (6) 

but the inner radius does not seem to scale. The latter fact was already mentioned in |||. The corresponding graph shows a 
sea-saw structure, i.e. short straight lines with jump discontinuities at certain values of the parameter p (see Fig. 5b). 




(a) (b) 


FIG. 5: rj = 0.19. (a) The locus of the outer apparent horizon is shown to satisfy the power law (6) with the slope 'Jdss = 0.109. The 
wiggles superimposed on the linear fit have the period 4.8 which agrees with the analytic prediction A/2y_D ss- (b) The locus of the 
inner horizon is shown not to scale. The jump discontinuities are periodic with the period 4.8. 

In order to understand this strange behavior we looked in more detail at the evolution of initial data fine-tuned to the 
location of these jumps. With the help of high resolution numerical methods we found the following remarkable structure: for 
a given family of initial data there is a sequence of discrete parameter values p\ > P 2 > P 3 -- > p n such that a solution with 
































5 


p € (jPmPn+i) approaches to the CSS\ solution n times, i.e. the solution comes close to the CSS 1 solution, turns away and 
returns n times before leading to black hole formation. Multiple approaches to the CSS\ solution were already noticed in 
|4) where they were called episodic CSS, however the corresponding fine-structure in the parameter space was not seen there. 
The sequence {p n } with n < 5 is listed in Table 1. 


n 

l 

2 

3 

4 

5 

oo 

Pn 

0.529001923689295 

0.528771570563995 

0.528769577618968 

0.528769560376615 

0.528769560227152 

0.528769560226138 


TABLE I: The first five critical parameter values p n for the coupling constant rj = 0.19. 


Because of numerical limitations we were not able to resolve higher p ni however the data shown in Table 1 seem to indicate 
that the sequence p n converges to p* as n tends to infinity. Actually, we find that the two consecutive parameters p n satisfy 
the scaling law 


Pn -P* _ ( A 

Pn+1 ~ P* ~ 6XP V2y dss) 


( 7 ) 


Now we return to the problem of scaling of the inner radius ry n . For p = p n + e, i.e. for p just above one of the p n 's we see a 
clear CSS scaling (see Fig. 6a) 


Tin ~ (p - Pn) 1CSS ■ (8) 

For p = p n — e the solution displays a kind of pseudo-dispersion after its last CSS-episode. This pseudo-dispersion manifests 
itself as follows: after leaving the CSS solution, the maximum of the function F decreases, the inner minimum of A disappears 
and later a spike develops which leads to the formation of an apparent horizon at r ou t er • In this range of p we see the subcritical 
CSS scaling (see Fig. 6b) 


max \d r F(t, 0)| ~ (p n — p) lcss , (9) 

however the masses of black holes formed in such evolutions are ’’large” and do not scale. We remark that since the solutions 
on both sides of p n form black holes, the bisection which gives critical parameter values p n has to be performed in a sense ”by 
hand”. 




(a) (b) 

FIG. 6: p = 0.19. (a) The supercritical (8) and (b) subcritical (9) scalings around p n for n = 2 (we get the same picture for each n). 
The slopes of the linear fits are equal to ±0.195 which agrees with the analytic prediction ±1/Acss where A css(v = 0.19) = 5.1 was 
obtained independently from the linear perturbation theory by Lechner §f]. 
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Interpretation of numerical results 

The results presented above confirm and extend the findings of Lechner et al. @]. Probing the bifurcation point r] c with 
higher accuracy we improved the evidence that A diverges as p tends to the critical value p c = 0.17 from above, which in turn 
confirms the picture that the DSS cycle merges with the CSS solution at the critical coupling constant p c . A natural question 
is: what is the meaning of the series of critical parameter values p n within this picture? 

We conjecture that our system shows a so called Shil’nikov bifurcation [12]. In his classification of loop bifurcations for three 
dimensional systems, Shil’nikov considered a system with a saddle point together with a homoclinic orbit which bifurcates for 
some value of a parameter. Assuming that the eigenvalues of the saddle point are real and satisfy the following conditions: 
Ai > 0 > A 2 > A 3 and Ai + A 2 >0 (plus some less important technical conditions), Shil’nikov showed that a saddle limit cycle 
bifurcates and the phase space picture looks qualitatively as in Fig. 7a. 




(a) 


(b) 


FIG. 7: (a) Shil’nikov bifurcation, (b) The conjectured phase space picture. 

Of course, our system is infinite dimensional and the Shil’nikov theorem cannot be applied directly. Nevertheless, it is 
expected that a similar picture to Fig. 7a will be valid for higher dimensional systems as long as only a few largest eigenvalues 
of the perturbation matter. Recently, Donninger ) 1 .‘li | has studied linear perturbations around the CSS solution and found 
that for coupling constants around the critical value p c the first three largest eigenvalues do in fact satisfy the above stated 
Shil’nikov conditions. Combining this property with the fact that the bifurcating DSS solution is a saddle limit cycle, we 
conjecture that the (one-dimensional) unstable manifold of the DSS solution lies on the stable manifold of the neighboring 
CSS solution. This is sketched in Fig. 7b. More precisely, the DSS unstable manifold winds around the limit cycle (infinitely 
many times) and eventually runs into the CSS saddle. Suppose that a curve of initial data intersects this spiral manifold at 
values p n , with limp n = p*, where p* corresponds to the intersection with the limit cycle. Then, the dynamical behavior will 
have exactly the form we observed above: for p equal to one of the p n ' s, the solution spirals n times around the limit cycle 
each time coming closer to the CSS-saddle before hitting it. For p = p n ± e, the behavior is similar, except that the solution 
does not hit the CSS-saddle but escapes along its unstable manifold. This is the reason why one observes the CSS scaling 
around p n with an exponent related to the unstable eigenvalue of the CSS solution. Note that the scaling law (7) follows 
immediately from the picture shown in Fig. 7b because during one cycle of evolution the distance from the DSS limit cycle 
increases by the factor e AossA (and 7 dss = 1/A dss)- 
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